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Abstract. Thermal photons are thought to be the ideal probe to measure the 
temperature of the quark-gluon plasma created in heavy ion collisions. PHENIX 
has measured direct photons with px < 5 GeV/c via their internal conversions 
into e+e~ pairs in Au + Au collisions at ^/sjvjv = 200 GeV and has now provided 
a baseline measurement from p + p data. 
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Ph i 1. Introduction 

»vj ! Direct photons are produced during all stages of heavy ion collisions at the Relativistic 

^ ' Heavy Ion Collider (RHIC). Because they do not interact strongly, they escape the 

0^ ! medium unaffected by final state interactions and provide a promising signature of 

^^ ' the earliest and hottest stage of the quark-gluon plasma (QGP) [T]. The main sources 

of direct photons from a QGP are quark-gluon Compton scattering {qg -^ jq) , quark- 
antiquark annihilation {qq —^ "fg) and bremsstrahlung involving thermalized and 

^^ I incoming partons [2 . At RHIC energies thermal photons are predicted to be the 

dominant source of direct photons in 1 < px < 3 GeV/c P . 

f~^ I The formation of a hot and dense medium has been established by experimental 

results at RHIC [3]. The measurement with the Electromagnetic Calorimeter 
(EMCal), based on a statistical subtraction of the expected background form hadronic 
sources, remains limited at low px due to large systematic uncertainties HUHj. 



H I The uncertainties in the knowledge of the decay background can be reduced by 

avoiding its explicit measurement. A tagging method has been developed in which 
a very pure photon sample is selected with strict photon identification cuts [5]. In 
this sample those photons are tagged, which combined with other clusters in the 
EMC (selected with less stringent photon identification cuts) can be identified via 
their invariant mass as the result of a 7r° decay. In the ratio of pure photons and 
tagged photons the reconstruction efficiency of the pure photons cancel and with it 
the associated systematic uncertainty. 

Furthermore, to circumvent the limitations due to the energy resolution at 
low photon energies, the excellent capabilities of the PHENIX detector to measure 
electrons can be utilized by measuring photons via their external conversion into e^e~ 
pairs [7]. These conversion pairs can be used as an alternative clean photon sample 
in the method described above. 
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2. Internal Conversions 



To overcome the statistical and systematic limitations an alternative approach has 
been pursued [8] and is outlined in the following. Any process that produces a real 
photon can also create a virtual photon which converts internally into an e~^e~ pair. 
Therefore, e^e~ pairs are also produced through internal conversions of virtual direct 



photons, e. g. via the qg Compton scattering {qg - 
between photon production and the associated e+e~ 
as: 



^ 1*9 ^ 6 6 9)- The relation 
pair production can be written 
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with TO being the mass of the e^e pair and ra^ ■ 
The process dependent factor S goes to 1 as m 



511 keV/c^ the mass of the electron. 
> or TO ^ Pt- For tt" and r/ decays, 
S is given by 5 = \F{m'^)\'^{l — m'^ /ml)^ [9] where mh is the hadron mass and F{m?) 
the form factor. For e+e~ pair masses approaching rah, the factor S goes to zero. 
While the measurement of real thermal photons suffers from a large background of 
hadron decays, measuring virtual photons allows to select a mass range, in which ir^ 
decays are suppressed due to the cut-off in S. 

The following analysis is based on 800M minimum bias events collected during 
the Au + Au run at ^/snn = 200 GeV in 2004 and 2.25 pb^^ of single electron 
triggered data recorded during the p + p run at -/s = 200 GeV in 2005. Electrons 
have been identified using information of the Ring Imaging Cherenkov detector (RICH) 
and a matching of momentum and the energy deposited in the EMC. Details of the 
analyses can be found in [10[ 111]. The invariant mass spectrum of e+e" pairs after 
background subtraction is shown in Fig. [1] for p + p (left) and Au + Au (right) for 
various e'^e~ pair-px bins above 1 GeV/c. The data are corrected for the electron 
identification efficiency based on a full Monte Carlo simulation, but are not corrected 
for the geometrical acceptance [TUl [TT] . These data are then compared to a cocktail 
of hadron decays from a fast Monte Carlo simulation whose input is based on the 
measured hadron spectra [101 E] ■ The cocktail is normalized to the data in the mass 
range 0-30 MeV/c^. While in p + p the agreement of the cocktail to the data is 
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Figure 1. Left: The e+e~ pair invariant mass distributions in (a) p + p and 
(b) minimum bias Au + Au collisions compared to a cocktail of hadronic sources. 
Right: Mass distribution of e+e~ pairs for Au + Au minimum bias events for 
1.0 < pt < 1-5 GeV/c with a fit as described in the text. 
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remarkable and a small excess is only observed at very high pt and large mass, the 
Au + Au data show a large excess over the full px range of 1-5 GeV/c. 

The observed excess above the hadronic cocktail can be analyzed under the 
assumption that it is solely due to internal conversions of direct virtual photons. 
While at zero mass e+e~ pairs from hadron decays have the same shape in mass as 
internal conversions of direct photons, the suppression due to S when approaching the 
hadron mass changes the shape of e+e^ pairs from hadron decays. Therefore, one can 
fit, after an initial normalization to the mass range 0-30MeV/c^, the two expected 
shapes for e"'"e~ pairs from hadronic decays /cocktail (as shown in Fig. [1]) and from 
direct photons /direct in the mass range 80~300MeV/c^ in which the tt" contribution is 
severely suppressed. The only fit parameter is the relative fraction of direct photons r: 

f{m) = (1 - r)/cocktaii("i) + rfdircctim). (2) 

As an example, the fit result is shown in Fig. [T]for e+e" pairs with 1.0 < px < 
1.5 GeV/c in minimum bias Au + Au collisions. The quality of the fit result (e. g., 
X^/NDF = 13.8/10 for the lowest pT bin) indeed justifies the assumption that the 
observed excess is due to internal conversion of virtual photons. A fit with the 
mass shape of e~^e~ pairs from 77 Dalitz decays would lead to a two times larger 
r] cross section than measured [12] and a significantly worse x^ /NDF = 21.1/10. 
Little contribution from other sources is expected to this kinematic region, as it is 
limited to m < 27tt,^. PHENIX has reported a large dielectron enhancement for 
150 < m < 750 MeV/c^ in [11], but the dominant fraction of the enhancement is 
concentrated at low pt [13] ■ 

3. Results 

The result is shown as 1 + r w 7inci/7hadr in the left of Fig. [2] for p + p collisions. It 
is compared to a NLO pQCD calculation of direct photons [H]. While the fraction 
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Figure 2. Direct photon excess, (a) The fraction of the direct photon component 
as a function of px in p + p. (b) Au + Au (min. bias) compared to other 
measurements of direct photons as described in the legend. The curves are from 
a NLO pQCD calculation [14]. 
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of direct photons measured in p + p agrees well with the pQCD calculation, a clear 
excess is observed in minimum bias Au + Au collisions shown in the right of Fig. [2 
The Au + Au result is also compared to the other measurements of direct photons 
at low pt; the conventional statistical subtraction of hadronic decay photons |15j . 
the tagging method [6] and the external conversion analysis [7]. While they all agree 
within their uncertainties, the improvement in statistical and systematic uncertainties 
is quite significant. 

The fraction of direct photons r can be converted into a direct photon spectrum 
by multiplying with the inclusive photon spectrum. The inclusive photon yield 
is determined for each px bin by dA^;"'^' = dN^f"" x dN^°''^*''^ / dN^°''^^'''\ where 
^jydata g^^^ ^^cocktaii g^j.g ^^^ yields of e+e~ pairs in m < 30 MeV/c^ for data 



and cocktail, respectively, and dNl 



cocktail 



is the yield of photons from the cocktail. 



The resulting direct photon spectra for p -\- p and Au + Au collisions are shown 
in Fig. [3l The p + p data can be well described with a modified powerlaw function: 
App{l +p^/6)~". The direct photons measured in Au + Au are above the binary scaled 
fit to the p + p result. An exponential fit plus the binary scaled modified power law 
(Ae^P-^^'^+TAA^App{l+p'^/b)'") reveals an inverse slope of T = 220±23±18MeVfor 
central collisions, suggesting an excess possibly of thermal origin. The inverse slope of 
the direct photon spectrum can be related to the intitial temparature of the QGP [16] . 
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Figure 3. Invariant cross 

section (p + p) and invari- 
ant yield (Au + Au) of di- 
rect photons. The p + p re- 
sult is compared to a NLO 
pQCD calculation. Solid 
points are from this analy- 
sis, open points from l4l[5l. 
The dashed lines show a 
Taa scaled, modified pow- 
erlaw fit of p-|-p. The black 
curves show an exponential 
plus the Taa scaled p + p 
fit. 
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